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Abstract 

We study the contribution of Galactic sources to the flux of astrophysical neutrinos 
recently observed by the IceCube Collaboration. We show that in the simplest model of 
homogeneous and isotropic cosmic ray diffusion in the Milky Way the Galactic diffuse 
neutrino emission consistent with 7 -ray (Ferrm-LAT) and cosmic ray data (KASCADE, 
KASCADE-Grande and CREAM) is expected to account for only 4%— 8 % of the Ice- 
Cube flux above 60 TeV. Direct neutrino emission from cosmic ray-gas (pp) interactions 
in the sources would require an unusually large average opacity above 0.01. On the 
other hand, we find that the IceCube events already probe Galactic neutrino scenar¬ 
ios via the distribution of event arrival directions. Based on the latter, we show that 
most Galactic scenarios can only have a limited contribution to the astrophysical signal: 
diffuse Galactic emission (< 50%), quasi-diffuse emission of neutrino sources (< 65%), 
extended diffuse emission from the Fermi Bubbles (< 25%) or unidentified TeV 7 -ray 
sources (< 25%). The arguments discussed here leave, at present, dark matter decay 
unconstrained. 



1 Introduction 


Galactic cosmic rays (CRs) interact with gas and radiation during the acceleration in their sources 
or during diffusion in the Galactic medium. Neutral pions produced in this interaction promptly 
decay into pairs of 7 -rays and are one component of the observed 7 -ray emission of the Galaxy. 
A significant component of the Galactic diffuse 7 -ray emission in the 100 MeV to 100 GeV range 
observed with EGRET [1| and Fermi-LAT [2] can be attributed to CR interactions with gas and 
is expected to become increasingly important towards higher energies. Direct evidence of hadronic 
interactions in Galactic sources could be inferred via the observation of low-energy cutoffs in 7 -ray 
emission spectra related to the production threshold of neutral pions [3]. 

The production of pions in the same CR interactions and the subsequent decay via 7 r + —> 
H + Vn followed by // + —» e + L' e v /J (and the charge-conjugate process) predicts a flux of high-energy 
neutrinos. The mean energy of a secondary neutrino produced in the production and decay of 
charged pions is typically of the order of 5% of the initial CR nucleon energy. Galactic CRs are 
believed to dominate the observed CR spectrum at least up to the CR knee, a spectral break at 
about 3 — 4 PeV. Therefore, for proton-dominated CR spectra the Galactic emission of 7 -rays and 
neutrinos is expected to extend to at least an energy of a few 100 TeV and possibly even higher. 

Interestingly, the IceCube Collaboration has recently identified a flux of high-energy neutrinos 
in the TeV-PeV energy range [4,5]. The analysis is based on high-energy starting events (HESE), 
i.e. events with neutrino interaction vertex inside the detector, as opposed to signals ranging into 
the detector. Figure 1 shows the arrival direction of a subset of HESE events with deposited 
energy of E ^ ep > 60 TeV in Galactic coordinates. This high-energy sub-sample shows a significant 
contribution of astrophysical neutrinos at the 5.7cr level. The best-fit isotropic F ” 2 power-law flux 
is at the level of E^dN/dE u = (0.95 ± 0.3) x 1CP 8 GeV cm -2 s ^ 1 sr -1 per flavor. 

The topologies of the HESE events are classified in terms of tracks and cascades, depending on 
whether the neutrino interaction produced a muon track inside the detector or just a nearly spherical 
emission pattern at its interaction vertex, respectively. In Fig. 1 these different event topologies are 
shown as diamonds and filled circles, respectively. The area of the symbols scales with the deposited 
energy of the event. The most energetic events are three PeV cascades (14, 20 & 35). At these high 
energies, absorption of neutrinos in the Northern Hemisphere due to charged and neutral current 
interactions as they pass through the Earth becomes relevant. As an illustration, the red-shaded 
area shows the minimal Earth absorption of neutrinos in the sample in log 10 -steps of 0.1. 

The IceCube Collaboration has analyzed a possible association of HESE events with known 
Galactic point-sources and a diffuse neutrino emission along the Galactic Plane; Galactic emission 
could not be significantly established [5]. On the other hand, an isotropic distribution of event 
arrival directions is expected for an extragalactic source population. As an illustration, the two 
plots of Fig. 2 show the HESE event distribution corresponding to events already shown in Fig. 1, 
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Figure 1: Mollweide projection of the arrival directions of IceCube events from Ref. [5] with 
deposited energies E^ep > 60 TeV in Galactic coordinates. The events are classified as tracks 
(diamonds) and cascades (filled circles). We use the same event numbers as in Ref. [5]. The relative 
detected energies of the events are indicated by the areas of the symbols. The thin lines around 
the arrival direction of the cascade events indicate the systematic uncertainty of the reconstruction. 
The red shaded region shows the minimal (E u = 60 TeV) absorption of the neutrino flux due to 
scattering in the Earth in log 10 -steps of 0.1. 




Figure 2: The Galactic latitude (left) and longitude (right) distribution of the HESE 3yr events 
with Efe P > 60 TeV. The blue solid lines show projections of the sum of probability distribution 
of event arrival directions normalized to an isotropic distribution. The red dashed and green dot- 
dashed lines show the expected distribution of an isotropic flux of 60 TeV and 1 PeV neutrinos, 
respectively, taking into account Earth absorption effects. The vertical positions of diamonds and 
circles indicate the deposited energy (right axis) of tracks and cascades, respectively. The dotted 
vertical lines indicate the central fit of the event location. 
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but now projected onto Galactic latitude (left) and longitude (right). In this plot we indicate 
the deposited energy via the height of the symbols, corresponding to the scale of the right axis. A 
strong emission from e.g. the Galactic Plane could be visible via an event cluster in the center of the 
latitude plot. However, the presently limited event statistics and the systematic uncertainties of the 
neutrino arrival directions make an identification of this Galactic structure less straightforward. In 
other words, the absence of a strong correlation of neutrino events with Galactic neutrino emission 
templates does not necessarily imply that a Galactic origin is already ruled out. 

Indeed, some authors have suggested that a significant contribution of this flux could have a 
Galactic origin (see Ref. [ 6 ] for a recent review). Possible sources are unidentified Galactic PeV 
sources [7,8] or microquasars [9], pulsar wind nebulae [10], extended Galactic structures like the 
Fermi Bubbles [11—14], the Galactic Halo [15] or Sagittarius A* [16]. A possible association with a 
hard diffuse Galactic 7 -ray emission [17,18] has also been investigated as well as a contribution of 
neutrino emission from decaying heavy dark matter [19-25]. A guaranteed flux of Galactic neutrinos 
is provided by the diffuse emission of Galactic CRs [12,26-33]. 

In this study we will discuss the expected high-energy emission of Galactic neutrinos with a 
special emphasize on its morphology. We provide estimates of the diffuse neutrino emission of the 
Galaxy in the TeV to PeV energy range utilizing the GALPROP propagation code [34], Our results are 
consistent with the hadronic sub-TeV 7 -ray emission inferred by the Fermi Collaboration [2], We 
will estimate the uncertainties on the overall flux from uncertainties of the CR spectrum in the knee 
region. The inferred emission rate density of Galactic CR sources can then be utilized to estimate 
the direct emission of 7 -rays and neutrinos, which can be parametrized by the opacity of the sources 
with respect to CR-gas interactions. We then try to estimate via an unbinned maximum likelihood 
test the sensitivity and discovery potential of the IceCube observation for Galactic emission. 

Our paper is organized as follows. We first calculate the Galactic diffuse neutrino flux in sec¬ 
tion 2, based on the GALPROP propagation code and with both power-law and exponential function 
fit to the CR spectra. In section 3, we discuss the fluxes and morphologies of point-like and extended 
sources. We then perform an anisotropy test of Galactic emission and calculate the sensitivity and 
90% upper limit on the Galactic fraction for different neutrino source candidates in section 4. Fi¬ 
nally, we conclude in section 5. Supplementary material for our study is provided in the appendices. 
In Appendix A, we compare our Galactic 7 -ray distributions with results of the Fermi Collaboration. 
We discuss the dark matter decay spectra in Appendix B, the background and signal samplings in 
Appendix C and the likelihood function in Appendix D. In Appendix E we discuss the extragalactic 
neutrino emission associated with other galaxies similar to the Milky Way. 
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2 Galactic Diffuse Neutrino Flux 


A guaranteed component of the diffuse 7 -ray emission of the Galaxy is due to CR interactions with 
the interstellar medium. Hadronic interactions with gas such as pp —> pp + 7 r° produce neutral 
pions that decay into two 7 -rays. The associated process pp —> pn + 7 r + can generate charged pions, 
which then decay into neutrinos. Therefore, the observation of diffuse Galactic 7 -ray emission in 
combination with the local measurement of CRs determines the flux and morphology of diffuse 
Galactic neutrino emission [12,26-33]. 

To calculate the secondary diffuse neutrino flux with an energy above 1 TeV, it is important 
to understand the primary CR flux, especially at energies around the CR knee at a few PeV. This 
feature could be related to the maximum rigidity achievable in Galactic accelerators or a change 
in the rigidity-scaling of diffusion. In both cases we expect a strong dependence on the chemical 
composition of Galactic CRs. In this study, we will adopt two Galactic CR models. The first one 
assumes a CR composition of proton and Helium following a broken-power-law spectra that we fit 
to the CREAM [35], KASCADE [36] and KASCADE-Grande [37] data. Our best-fit CR spectrum 
follows 
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and E^ = 8.2 PeV. In the left panel of Fig. 3, one can see that our power-law fit in good agreement 
with to the observed CR data up to the knee region. Note, that CR diffusion in the Galactic medium 
softens the observable spectra. As a consistency check of our implementation of the CR injection 
spectrum Eq. (1) into GALPROP one can compare the energy dependence of the resulting CR flux to 
analytic expectations. The difference between the observed CR power E~@ and the injected power 
E~ a is expected to follow the relation (3 = a + 6, where 5 is the rigidity-dependence of the diffusion 
coefficient. In this study, we assume diffusion of CRs in a Kolmogorov-like magnetic turbulence that 
predicts 5 = 1/3. The numerical prediction of GALPROP shown in the plot reproduces this expecting 
rigidity-scaling well. 

As a second CR model we adopt the best-fit result of a study by Gaisser et al. [38] modeling 
the CR spectrum as a superposition of three populations with individual contributions of different 
mass groups that have a common exponential rigidity cutoff [39]. As in the first CR model, the 
proton and helium components still provide the dominant contribution to the neutrino flux. The 
CR spectrum is given by a parametric form in Ref. [38] [Table II and Eq.(3) therein] and is shown 
in the right panel of Fig. 3. The parametrized spectra are given by the sum over three populations, 
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Figure 3: Left panel: a broken-power-law fit to the CR spectrum including only proton and 
helium contributions. Right panel: a more sophisticated fit with exponential functions from 
Ref. [38], where the iron flux is important around the knee. 

with normalization a p = (0.34, 0.020, 0.00032) GeV m -2 s -1 sr -1 , cutoff energies E* = (4, 30, 
2000) PeV and spectral indices = (2.66,2.4,2.4) for the three proton populations and an e = 
(0.35, 0.015, 0.00025) GeVm -2 s -1 sr -1 , E^ c = 2 E* and Tne = (2.58, 2.4, 2.4) for helium. The 
two CR models have different powers below the knee and the power-law model provides a harder 
spectrum, i.e. larger CR fluxes. In the low energy part below 100 TeV, neutrino fluxes follow a 
simple power law behavior and have a power of — 2.54(—2.69) close to the proton CR fluxes in 
Eq. (1) [Eq. (2)]. 

In the left panel of Fig. 4, we show our predictions for the all-sky-averaged diffuse neutrino fluxes 
for the two different CR models. We also show the best-fit power-law flux observed by IceCube [5], 
E^dN/dEv = 4.5 x 10 -8 (I%/100 TeV) -0 ' 3 GeVcm _2 s _1 sr -1 summed over neutrino flavors. The 
diffuse neutrinos account for 10% (5%) for the broken-power-law (global) fit with the sky-averaged 
flux and E v i v ^ 30 TeV (see Refs. [40-44] for other spectral analysis). For E u / V ^ 60 TeV, the diffuse 
neutrinos account for a slightly smaller fraction of 8% (4%) for the broken-power-law (global) fit. 
Compared to the atmospherical muon neutrino background measured by IceCube [45] (see also 
[46-48] for theoretical calculations), the Galactic diffuse neutrino flux dominates over the Galactic 
diffuse spectrum below PeV. However, this comparison is based on neutrino fluxes averaged over 
the whole sky. As we know from the density distribution indicated in Fig. 11 of Appendix A, 
the Galactic diffuse neutrinos mainly come from the Galactic Plane. In the right panel of Fig. 4, 
we show the comparison for the neutrino differential fluxes around the Galactic Plane region with 
\b\ ^ 7.5° [45]. The Galactic diffuse neutrinos can have a flux over the atmospheric one for energy 
above ~ 300 TeV. We also note that the IceCube HESE analysis (see Fig. 1) has reported a reduced 
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Figure 4: Left panel: the differential neutrino and 7 -ray spectra averaged over all the sky, 
based on the GALPROP s S Z 4 R 20 1 150^5 diffuse model. Also shown is the IceCube best-fit and 
all-flavor spectrum flux: E^dN/dE v = 4.5 x 10 _8 (£ , J/ /100 TeV)~ a 3 GeVcm - 2 s - 1 sr -1 [5]. For the 
7 -ray spectrum, absorbing effects become significant only above 10 TeV due to the CMB absorbing 
effects [57]. We also show the IceCube measured atmospherical muon neutrino background flux in 
the green and cross points [45]. Right panel: the same as the left one but for the Galactic Plane 
region with \b\ ^ 7.5°. 

atmospherical neutrino background and decmonstrated an enhanced ability to measure Galactic 
diffuse neutrinos [49]. 

Before we move on to the next section, a few remarks are in order. Our analysis of Galactic 
diffuse gamma-ray and neutrino emission with the GALPROP code follows the standard approximation 
of homogenous and isotropic CR diffusion in the Milky Way. It also assumes that the observed CR 
spectrum can be approximated as a steady-state solution of a smooth distribution of continuous 
sources. All these assumptions can only be considered as a first-order approximation. Deviations 
from these assumptions can decrease or enhance the diffuse emission. For instance, it is well known 
that diffusion along the regular component of the Galactic magnetic field is stronger compared 
to orthogonal directions, e.g. [50]. This can lead to local variations of the CR flux normalization 
compared to the Galactic average [51]. A similar effect can be observed for models accounting for 
non-azimuthally symmetric source distributions including the Galactic spiral arms. This has been 
studied with three-dimensional propagation codes like DRAGON [52,53] and PICARD [54,55]. It was also 
shown that non-homogenous Galactic diffusion implemented in the propagation code DRAGON [52] 
can enhance the hadronic gamma-ray and neutrino emissions in the multi-TeV region [56]. Finally, 
it was argued that time-dependent and local CR injection episodes could also result in locally 
observed CR spectra that are softer than the Galactic average [17]. This harder average emission 
would increase the diffuse emission in the multi-TeV region after normalization to Fermi data. 
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3 Neutrino Fluxes from Point-like and Extended Sources 


Point-like sources in the Milky Way may also contribute to the observed neutrino flux. Since there 
is presently no evidence for a correlation of the IceCube HESE events with Galactic point-like 
sources [58] we do not expect a strong contribution of individual Galactic sources, but rather look 
for hints in the data for their combined (stacked) emission. The most likely candidates of this 
neutrino emission are the sources that contribute the bulk of Galactic CRs. Candidate sources 
include supernova remnants (SNRs) and pulsar wind nebulae (PWNe). 

To first order, the observed density ticr of CRs can be related to the time and volume-averaged 
CR emission rate of the Galaxy (QcrX-^ 0 via the energy-dependent diffuse escape time T esc (E ) as 
ncR — T esc (QcR){E) ( see e -9- Ref. [59]). The escape time from the Galactic Plane with half height 
H can be estimated as r esc — H 2 /D(E) for a diffusion coefficient D(E). The initial spectrum of 
CRs is therefore softened by the energy dependence of D(E ). 

The fit to the observed CR data via GALPROP presented in the previous section can also be 
utilized to determine the in situ contribution of the sources of Galactic CRs more precisely. From 
the GALPROP fit we know the time-averaged differential source emission rate of Galactic CR sources 
Q cR(r, z, E) for a given source distribution and propagation model. The difference to the diffuse 
neutrino emission of Galactic CRs is that here we assume that secondary neutrinos are directly 
produced at the source location via interaction with gas in the source environment, e.g. Refs. [60,61]. 
This is quantified via the optical thickness T pp . So, one can directly relate the source emission rate 
of neutrinos (per flavor) and CRs as 


1 


ElQ v {r, z, E v ) * - (1 - e~ KineTpp ) E 2 Q C n(r, z, E p ). (3) 

with inelasticity Ki ne ^ 0.5. The overall factor of 1/6 is a rough estimation of the total neutrino 
energy over the total proton energy accounting for the per flavor emission (1/3), the total neutrino 
energy fraction in the charge pion decay (3/4) and for the charged pion fraction in pp collisions (2/3). 
Furthermore, for the individual neutrino energy, one can use the approximate relation E u ~ E p /20 
with a factor of 1/4 for the neutrino (each flavor) energy in the pion decay and a factor of 1/5 for 
the average energy of the pion produced in pp interactions. 

For a source with dynamical time-scale fdyn and average gas density n gas the opacity is given 
as T pp ~ ctdyn^gasCpp for an inelastic scattering cross section a pp ~ 34 + 1.88L + 0.25L 2 mb for 
L = ln(Sp/TeV) [62,63]. We estimate the dynamical time-scale as 10 4 yrs corresponding to the 
beginning of the snowplow phase, where the acceleration efficiency in typical supernova remnants 
are expected to cease [64], We therefore estimate the optical thickness, following Refs. [60,61] as 
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For a given class of sources, the approximate Galactic source distribution is typically expressed in 
cylindrical Galacto-centric coordinates p(r, z), which is related to spherical Helio-centric coordinates 
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(s, £, b ) as r 2 = s 2 cos 2 b + Rq — 2 sRq cos l cos b and z = s sin b. The distance from the Sun to the 
Galactic Center is Rq ~ 8.5 kpc. To keep track of the source distribution in the skyrnap, it is 
convenient to introduce the normalized “ J”-factor after integrating along the line of sight 

J(b,£) = ^dsp[r(s,£,b),z(s,£,b)] . (5) 

with the total number of sources N s = 2 r irS t dzS l drr p(r,z). The total J-factor is defined to be 
J = ^dLlJ(b,£) = \ dftdfcos b J(b, £). We can then express the all-sky-averaged neutrino diffuse 
neutrino flux via the total neutrino spectral emission rate Q u as 


R u 4>u 
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( 6 ) 


3.1 SNRs and PWNe 

In the following, we study two CR source densities that are expected to approximate the distributions 
of SNRs and PWNe in our Milky Way. For SNRs we use the parametrized distribution from Case 
et al. in [65], while for the PWNe, we take the pulsar distribution from Lorimer et al. in [ 66 ]. Both 
can be described by the following function 

p(r ’ 2) = (jfe) exp (~ pL fl?) exp (“x) ’ (7) 

with a = 2, [d = 3.53 and h = 0.181 kpc (Case, SNR) and a = 1.93, /3 = 5.06 and h = 0.181 kpc 
(Lorimer, pulsars). Both models give J ~ (16.2 kpc)^ 2 . 

In Fig. 5 we show the estimated all-sky-averaged neutrino flux from CR sources in our Galaxy 
assuming a distribution following that of SNRs. We assume an average optical thickness of the 
sources of r pp ~ 3 x 10~ 4 following the SNR benchmark values shown in Eq. (4). Note that the 
flux level is proportional to the source opacity as long as r pp « 1 and saturates for r pp » 1. To 
have the flux not overshoot the observed IceCube flux, the averaged opacity is required to be 
< \t pp ) < 0.01. The overall normalization is determined from the GALPR0P fit to CR data assuming 
two different fitted CR spectra in Eq. (1) and Eq. (2) resulting in a total Galactic CR power of 
about (3 — 5) x 10 10 erg/s. This is consistent with a Galactic supernova rate of 3 yr _1 if a few 
percent of the kinetic energy of supernova ejecta of typically 10 51 erg is converted into CRs. Since 
the spectrum is harder than the diffuse emission by the diffusion parameter 5 = 1/3, the source 
spectrum can dominate at high energies unless the opacity of the sources becomes too small. 

3.2 Unidentified TeV Sources 

Other than the cumulative contribution of SNRs and PWNe, we also consider the contribution from 
unidentified TeV (UnID TeV) 7 -ray sources [7], which may originate in hypernova (superluminous 
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Figure 5: The data points are the 7 -ray flux for one representative UnID TeV source: MGRO 
J1908+06 [67,68], after multiplied by a factor of 1/47T. The solid lines are calculated neutrino fluxes 
directly from the super novae remnant sources. 


supernova) remnants of our Galaxy [69] and contribute to high energy neutrinos from hadronic pp 
interactions. For instance (c/. data points in Fig. 5), one of the bright sources, MGRO J1908+06, 
observed first by the Milagro collaboration, has the differential 7 -ray flux of, E 2 dN/dE = 1.3 x 
10 _9 (£'/1 GeV )^ 0 20 GeVcm" 2 s _ 1 sr _1 (after multiplying the source flux by l/ 47 r), from a more 
precise measurement by the VERITAS collaboration [ 68 ] which is just one order of magnitude below 
the IceCube neutrino flux. Two of the UnID TeV sources studied in [7] have since been identified 
(HESS J1018-589 & HESS J1837-069) whereas one new unidentified source has been announced 
(VER J2019+368). Consequently, we have updated the UnID source catalogue from the tabulated 
list provided in Ref. [7]. 

3.3 Fermi Bubbles 

Beyond the point-like sources, there are also a few candidates of extended neutrino sources. The 
Fermi Bubbles [70] are two large structures extending to latitudes of about 55° above and below the 
Galactic Plane. The precise mechanism of the 7 -ray emission is unknown and, both, hadronic [71] 
and leptonic scenarios [72, 73] have been proposed in the literature as an explanation. A possible 
connection to the IceCube HESE flux has been pointed out in Refs. [11-14], The Fermi Collaboration 
has also performed a dedicated search to understand the energy spectrum and morphology of the 
Fermi Bubbles [74], For the hadronic origin, a power-law CR proton spectrum with an exponential 
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cutoff at about 15 TeV provides a better fit than a power-law spectrum without a cutoff [74], The 
corresponding cutoff in the neutrino emission would then be below TeV. However, the significance of 
this cutoff is weak and it is possible that the emission continues with a softer power towards higher 
energies and then contributes to the IceCube data above 10 TeV [74], Future 7 -ray observations 
with CTA [75,76] and HAWC [77,78] in the (0.1 — 100) TeV range will be able to identify the 
spectrum [12,14], 


3.4 Decaying Dark Matter 


Another more exotic scenario that has been proposed as a source of the IceCube flux is decaying 
dark matter with a dark matter mass of 0(5 PeV) and a lifetime of 0(1O 28 s) [19-21,24,25]. The 
contributions of Galactic emission from the Galactic dark matter halo and an extragalactic isotropic 
emission are of similar magnitude [20,25,79,80]. The relative intensity of the Galactic emission is 
again calculated from Eq. (5). For the dark matter density distribution in our galaxy, we use the 
isotropic Einasto profile [81] with 


Pdm(t) = Pq exp 





( 8 ) 


with pq = 0.4 GeV/cm 3 , r s = 20 kpc, tq = 8.5 kpc and (3 = 0.17. The total dark matter mass in 
the Milky Way is then approximately ~ 2.7 x 10 12 Mq. 

Note, that the contribution from extragalactic dark matter will be mostly isotropic, except for 
close-by galaxies and galaxy clusters [25]. This isotropic background and the very extended Galactic 
diffuse neutrino emission will make it difficult to identify the emission by anisotropy studies alone. 
However, the dark matter decay spectrum can be identified by spectral properties, like strong line- 
features and mass thresholds that are unexpected in astrophysical scenarios. Two representative 
dark matter decay spectra are shown in Fig. 12 of Appendix B. 


4 Anisotropy Test of Galactic Neutrino Emission 

In general, Galactic neutrino emission is expected to be anisotropic. The diffuse emission from 
CR propagation is required to follow the integrated column density along the line-of-sight. This 
emission is structured due to the higher gas concentration towards the Galactic Plane (GP) and 
Galactic Center, but it also has high-latitude fine-structure due to molecular gas clouds. Emission 
from CR-gas interactions in the vicinity of CR sources is also correlated with the Galactic Plane. 
In the following we will consider two source density distributions already introduced in the previous 
section, approximating supernova remnants (Case [65]) and pulsar wind nebulae (Lorimer [ 66 ]). On 
the other hand, a hypothetical signal from the decay of heavy dark matter (DM) is expected to 
follow the DM density distribution producing an extended emission around the Galactic Center. As 
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Galactic diffuse (logjo( M 'signal/ M 'iso)) SNR (loglO^signal/v^iso)) 







Figure 6: Six Mollweide projections of Galactic emission templates used in the statistical analysis. 
From left to right and top to bottom, they are Galactic diffuse (E v = 10 TeV), supernovae remnants 
(Case), pulsar wind nebulae (Lorimer), decaying dark matter, Fermi Bubbles and unidentified and 
dark TeV 7 -ray sources. The mesh indicates the equatorial coordinate system with right ascension 
a = 0° and declination <5 = 0° indicated as solid lines. Note that at the location of IceCube 
declination and zenith angles have the simple relation 9 zen = 7t/2 + 5. 
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a template we use the Einasto profile [81]. To test for extended diffuse neutrino emission from the 
Fermi Bubbles, we use the gamma ray emission template provided by the Fermi Collaboration [74], 
The collective emission from Galactic hypernovae [7] is already covered by a test of emission along 
the Galactic plane. In addition to that we also include the unidentified and dark TeV 7 -ray sources 
that have been summarized in Ref. [7] which we update with recent observations. In Fig. 6 , we 
show the geometric distributions of all six Galactic emission templates in the Galactic coordinate. 
For Galactic diffuse, supernovae remnants (Case), pulsar wind nebulae (Lorimer) and decaying dark 
matter templates, we show their flux density distributions with respect to an isotropic distribution. 

Correlations of the neutrino arrival direction with Galactic structures can be tested using an 
unbinned maximum likelihood test (see e.g. [82]). The IceCube Collaboration has tested emission 
along the Galactic Plane via a top-hat template, wgp(Q) = @( 6 qp ~~ |&|)/(47rsin&Gp)> where + 6 qp 
is the Galactic lattitude width of the emission region. A scan over 6 gp revealed a maximal excess 
over background at 6 gp = 7.5°, with post-trial significance of only 2.8% [5]. 

Here, we want to improve on previous tests in several ways. Firstly, the likelihood ratio in 
our statistical test accounts for the intrinsic anisotropy of atmospheric backgrounds and Galactic 
emission as well as the expected anisotropy due to Earth absorption and detector acceptance. This 
intrinsic dependence on zenith angle does not cancel in the likelihood ratio and is not accounted 
for in the IceCube analysis. The likelihood ratio with respect to the background hypothesis (a full 
isotropic astrophysical signal) can be written in the form 

£(/iso) = FT 

£(1) -H 

where the signal and background expectation values /a S j g and /^j gr , respectively, for the j’th event 
depend on the fraction /; so of the total number of expected events from the best-fit FA 2,3 spectrum 
that is due to an isotropic emission (see discussions in Appendices C and D). The complement 
/Gal = 1 — /iso is the event fraction associated with Galactic sources. The maximum log-likelihood 
(LLH) ratio 

p (10) 

with maximum point /i so is then defined as our test statistic (TS), which will be used to estimate 
the sensitivity of the IceCube Observatory for Galactic emission and the significance of emission in 
the HESE analysis and a classical up-going search. 

For the sensitivity and significance studies of Galactic emission in the HESE data and the 
classical up-going z/ /t data we generate event maps following the prescription outlined in Appendix C. 
The left plots in Fig. 7 show examples of samples for the full signal case (/j so = 0, top), a partial 
distribution (/i so = 0.5, middle) and the background case (/i so = 1, bottom). The title of the plots 
indicate the test statistics value A and the value /j so at maximum. 


^ Slg (/iso) + h- gr (/-) 

!C Sr (l) 


(9) 
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Samples for 3yr HESE (E^ep > 60 TeV) 


Samples for 3yr Up-going (E^ > 100 TeV) 



o/o : Galactic v \ <>/o : isotropic v \ <>Jo : atmospheric v \ o/o : atmospheric p. 


Figure 7: Event samples with decreasing Galactic diffuse fraction /; so of 0 (top), 0.5 (middle) and 1 
(bottom) for the HESE 3yr study with E ^ ep > 60 TeV (left panels) and a classical up-going v tJ search 
with En > 100 TeV in the same period (right panels). In the sample maps the simulated events 
are composed of Galactic neutrinos (black), isotropic neutrinos (green), conventional atmospheric 
neutrinos (blue) and atmospheric muons (magenta). Tracks are indicated as diamonds and cascades 
as filled circles. The angular uncertainty (sampled from data) is indicated as circles. The numbers 
in all maps follow the relative signal strength of the events for a Galactic diffuse origin starting from 
strong to weak. Each map indicate total number of events ntotj the maximum likelihood ratio A 
and the corresponding maximum fraction f[ so . 
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Figure 8: Test Statistic (TS) distributions of 5000 Monte Carlo samples for the HESE 3yr study. 
For each model the histograms show background (gray) and 100% Galactic emission (red). The 
vertical red dashed line shows the median TS of simulated Galactic emission and the vertical green 
solid line the observed value Ahese- The black points indicate the bin-wise contribution of the 
expected distribution following [5(A) + Xi(A)]/2 [83]. 
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Figure 9: Sensitivity (left panel) and 90% upper limit (right panel) on the Galactic fraction 
/Gal = 1 — /iso for six different Galactic emission templates. 

Fig. 8 shows the TS distribution of 5000 samples generated for the background case (/; so = 1) 
and a full contribution (/i so = 0). The test statistics of the HESE 3yr data are indicated as vertical 
green lines and summarized as the second column of Table 1. For all six emission models we find that 
the data are consistent with the background distributions, i.e. an isotropic arrival distribution of the 
astrophysical signal. This is consistent with previous findings of the IceCube Collaboration [5]. The 
highest excess above background is obtained with the PWN distribution with a chance probability 
of 9%. 

The medians of the TS distributions (vertical red dashed lines) are well separated from the back¬ 
ground distribution in most cases. This indicates that the non-observation of significant Galactic 
anisotropies in the HESE 3yr data (vertical green line) is constraining the contribution of most 
Galactic emission scenarios. The 90% sensitivity to a Galactic fraction is defined as the minimum 
isotropic fraction /i so where 90% of the signal samples have a TS larger than the median TS of the 
background distribution [84], For the likelihood defined via Eq. (9) the median for all six Galactic 
emission templates is 0. We estimate the sensitivity via a scan over f\ so with A/; so = 0.05 steps and 
show our results in Fig. 9 and the sixth column of Table 1. 

All scenarios considered in this paper have been suggested before to explain a part of the HESE 
data. Therefore our study does not correspond to a “blind” analysis and it is no surprise that 
all models show an up-fluctuation above background. The 90% upper limit from this analysis 
corresponds to the maximum Galactic fraction where 90% of the signal samples have a TS larger 
than the value observed in the data. We derive this upper limit via a scan over Galactic contributions 
/Gal = 1 — /iso and show the result in the right panel of Fig. 9 and the fifth column of Table 1. 
Interestingly, for the case of Galactic dark matter decay we can not obtain a (non-trivial) upper 
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Figure 10: Sensitivity of the search of up-going above muon energies of 100 TeV in three years 
(left) and the sensitivity of the HESE search after 10 years of observation (right). In the case of the 
classical up-going search we do not include the templates of Fermi Bubbles and unidentified or 
dark TeV 7 -ray sources, since the overlap with the Northern Hemisphere is small. 


limit on the Galactic fraction based on the HESE 3yr data and hence we leave the entry empty in 
Table 1. 

With our method we can also estimate IceCube sensitivities for the various Galactic scenarios 
with increased statistics and with alternative up-going muon neutrino searches. For the classical 
up-going ( 9 zen > 85°) + search we assume a total astrophysical event rate of 10 events per year 

and a total conventional atmospheric background rate of 10 events per year. This approximates the 
expected integrated event rate above muon energies of 100 TeV [85], comparable to the cut of 60 TeV 
in deposited energy used for the HESE search. The details of the expected signal and background 
distributions are described in Appendices C and D. The right plots in Fig. 7 show examples of 
samples for the full signal case (/i so = 0, top), a partial distribution ( f lso = 0.5, middle) and the 
background case (/j so = 1, bottom) over a period of three years. Note that this approach is not 
sensitive to contributions which dominate in the Southern Hemisphere and therefore we leave out 
the Fermi Bubbles and UnID sources in this study. 

The left plot in Fig. 10 shows the result of the sensitivity scan for the muon neutrino search 
and the last column of Table 1 summarizes the results. Note that the sensitivities are not much 
different compared to the HESE 3yr search for the four scenarios, although the muon search offers a 
better angular resolution of the order of 0.5°. The reason for this lies in the fact that the templates 
themselves correspond to extended emission such that the improved pointing precision does not 
matter that much. The right plot in Fig. 10 shows the result of the sensitivity scan for 10 years of 
HESE observations and the sensitivity level is also summarized in the seventh column of Table 1. 
The sensitivity compared to the three years data improves by about a factor of two in this period. 
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HESE 3yr observation 


sensitivity for / Ga l* 

template 

A 

p- value* 

/Gal* 

f90%* 

J Gal 

HESE 

3 yr 

HESE 

10 yr 

Northern 
^ 3yr 

Galactic diffuse v # 

0.74 

0.19 

0.19 

0.50 

0.30 

0.15 

0.25 

SNR [65] 

1.68 

0.10 

0.34 

0.65 

0.35 

0.20 

0.30 

PWN [66] 

1.77 

0.09 

0.30 

0.60 

0.30 

0.15 

0.25 

DM decay [81] 

1.48 

0.11 

0.46 

- 

0.60 

0.30 

0.85 

Fermi Bubbles [74] 

0.36 

0.27 

0.07 

0.25 

0.20 

0.10 

- 

UnID TeV [7] 

0.43 

0.25 

0.07 

0.25 

0.20 

0.10 

- 


* The emission template is using GALPROP. We estimate the systematic uncertainty of /cal from the diffusion model to be 
at the level of +10%. 

* The p -value is calculated from A assuming a background distribution [<5(A) + x|(A)]/2. 

* The Galactic fraction is defined as /cal = 1 — /iso- 

Table 1: Sensitivity and 90% C.L. lower limits of a Galactic fraction in the HESE data above 
60 TeV. The first two columns shows the TS and maximum point n s using the IceCube approach 
via Eq. (9). 

Note that our result is not a full replacement of an IceCube analysis. Several steps of this 
analysis can be improved, in particular the zenith and energy dependence of the events. We expect 
that a dedicated IceCube analysis will improve the sensitivity of the analysis by a factor of a few. 
In particular, for very high energy neutrinos the classical muon neutrino is also sensitive to emission 
in the Southern Hemisphere, although at a much lower level [58]. A strong Galactic contribution 
can also alter the best-fit value of the astrophysical contribution which requires a simultaneous fit 
in the first place. 

5 Discussion and Conclusions 

We have studied the contributions of extended Galactic TeV-PeV neutrino emission sources in 
relation to the IceCube observations. A guaranteed contribution to Galactic emission is from CR 
propagations and interactions in the Galactic medium. We have studied the corresponding diffuse 
emission of gamma-rays and neutrinos with the numerical cosmic ray propagation code GALPROP. In 
our calculations we have assumed that the locally observed CR flux corresponds to the steady-state 
solution of the diffusion-convection equation with a homogeneous and isotropic diffusion coefficient. 
We found that under these assumptions the expected Galactic diffuse neutrino emission that is 
consistent with 7-ray (Ferrm-LAT) and CR data (KASCADE, KASCADE-Grande and CREAM) 
can only contribute 4%-8% above 60 TeV. We also derived the neutrino emission from the Galactic 
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CR sources parametrized by the opacity of the sources to CR-gas interactions. We showed that the 
same CR model assumptions require an unusually large average opacity of (r pp ) ^ 0.01 to reach the 
flux level of IceCube integrated over the whole sky. 

However, a more realistic study of CR injection and diffusion can alter these results, for instance 
time-dependent and inhomogeneous sources as well as inhomogeneous and/or anisotropic diffusion 
in the Galaxy. For instance, it was argued that remnant CRs from a locally enhanced CR emission 
episode could soften the CR spectra observed today compared to the global Galactic average [17]. 
Anisotropic diffusion in our local environment would similarly introduce local fluctuations from the 
Galactic average leading to, both, enhancements or reductions of the expected all-sky hadronic 
gamma-ray and neutrino emission [51]. A similar effect has been discussed in the context of more 
realistic source distributions following the Galactic Arms [54,55]. The combined effect can reach to 
local enhancement and deficits of the CR spectrum at the level of 0.5 — 3.0. Note, however, that 
these studies do not account for the Fermi gamma-ray data to constrain their model parameters 
and the actual range might be somewhat smaller. Finally, inhomogeneous diffusion could enhance 
the local concentration of CRs [52], It was argued in Ref. [56] that this would also increase the 
multi-TeV hadronic emission of the steady-state solution in the inner Galaxy by a factor of two in 
better agreement with Fermi and Milagro data. 

An independent limit on Galactic neutrino emission that does not rely on the overall normal¬ 
ization scale of CR injection can be derived from the absence of correlations of neutrino arrival 
directions with Galactic emission templates. We studied the published IceCube data with respect 
to Galactic diffuse neutrino emission, emission from Galactic source populations like SNRs or PWNe, 
emission from the Fermi Bubbles or unidentified/dark TeV 7-ray sources and a dark matter halo. 
We did not find a strong statistical excess in the 3yr HESE data above background expectations. 
This conclusion is consistent with previous studies of the IceCube data. However, we showed that 
most of these scenarios are already constrained by the 3yr HESE data. The contribution at the 90% 
confidence level is limited to < 50% for diffuse Galactic emission, < 65% for quasi-diffuse emission 
of neutrino sources and < 25% for extended diffuse emission from the Fermi Bubbles or unidentified 
TeV 7-ray sources. Interestingly, the emission of PeV dark matter decay in our Galactic center is 
presently unconstrained by the data. 

We also estimated the sensitivity of IceCube for these emission scenarios with 10 years of HESE 
data and with 3 years of a classical up-going muon neutrino search. We estimate that the classical 
muon neutrino search in the same time period (three years) and energy range (muon energies above 
100 TeV) has a similar sensitivity to extended diffuse Galactic emission in the Northern Hemisphere. 
The sensitivity of the HESE data after ten years of observation is expected to increase by a factor 
of two, not accounting for an up-fluctuation of the test statistic in the first three years. 

Our analysis only relies on published IceCube data and can not be considered a full replacement 
of a dedicated IceCube study. We expect that our sensitivity estimates and upper limits are conser- 
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vative and can be improved for the same data set using additional event information. In particular, 
our maximum likelihood analysis does not account for the energy dependence of individual events 
and can be improved by IceCube: if the Galactic emission spectrum responsible for the 60 TeV 
to 3 PeV data shows strong spectral variations with respect to the isotropic emission, a statistical 
separation between Galactic and extragalactic components can be much stronger. In particular, 
this is expected for dark matter decay, since narrow spectral features ( e.g . line features) of c.lose-by 
Galactic dark matter are expected to be smoothed out for extragalactic contributions due to the 
redshift distribution. 

Finally, let us conclude with a few additional remarks. The data-driven limit on Galactic 
diffuse neutrino emission is based on a template derived with GALPROP, assuming an azimuthal 
distribution of CR sources and homogeneous isotropic diffusion. As noted earlier, the template does 
not depend on the absolute normalization, but it is expected to vary for models that predict relative 
enhancements of local CR distributions [17,51,55,56]. This introduces a systematic uncertainty for 
the limit on Galactic diffuse neutrino emission. However, due to the dominance of cascade events in 
the IceCube data, which have a poor angular resolution of more than 10 degrees, we don’t expect 
that the statistical analysis of the HESE data has a very strong dependence on the diffusion and 
source model. 

One possibility to estimate the effect of local CR enhancements by source and diffusion models 
is by a comparison of the statistical outcome for the SNR and PWN templates. The relative ratio 
of the corresponding azimuthal source distributions (normalized to the same number of sources) 
is comparable to the variation of 0.5 — 3.0 shown in detailed studies [55] for non-azimuthal source 
distributions. The resulting templates shown in Fig. 6 differ in the maximum-to-minimum ratio by 
a factor two. However, from Fig. 9 and Tab. 1 we can see that the different statistical results for 
the SNR and PWN emission templates are within +10 % and we take this as an estimate for the 
systematic error associated to the modeling of Galactic CR emission and diffusion. 

A strong neutrino emission from our own Milky Way also implies an extragalactic contribution 
from similar galaxies. In Appendix E we show that this extragalactic contribution is in general 
expected to be not as significant as the Galactic emission. The situation is different for the case of 
PeV dark matter decay: Here, the extragalactic flux is comparable to the Galactic flux. In addition, 
any strong spectral features from particle physics, that are possible for the local Galactic emission 
in this scenario (see Appendix B), are expected to be smoothed out due to integration over red-shift. 

Hadronic production of TeV-PeV neutrino sources will inevitably predict TeV-PeV 7-rays [12], 
which can be observed by large-scale 7-ray telescopes like HAWC [14,77,78] or CTA [75,76]. In par¬ 
ticular, the observation of PeV 7-rays would correspond to a smoking gun of Galactic contributions 
due to the small absorption length of about 10 kpc via e + e _ production in the CMB [12,86-88]. 
For pp interactions, the relation between the 7-ray and neutrino fluxes is shown in Fig. 4, which 
has small dependence on the primary CR index. For decaying dark matter models, a wide range of 
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dark matter models predict the associated 7 -ray flux above 10% of the observed IceCube neutrino 
flux as shown in Fig. 13 in Appendix B. This is within the reach of HAWC-100 or HAWC-300 with 
one year of data [77]. 

In our study we have simulated Galactic and extragalactic (isotropic) event distributions as¬ 
suming an equal neutrino flavor ratio with (/ e : /), : f T )^ « (1 : 1 : 1)0 which corresponds to an 
initial flavor combination ( 1 : 2 : 0 ) s (pion decay) at the astrophysical sources, after flavor oscilla¬ 
tions [89,90]. This is consistent with IceCube observations [5,43,91,92]. Future observations could 
be used to distinguish different explanations for the observed IceCube events [93-97], in particular 
( 0 : 1 : 0 ) s (muon-damped pion decay), ( 1 : 0 : 0 ) s (neutron decay) and ( 0 : 0 : l) s (flavor-dependent 
dark matter decay). This fact makes our morphology studies of prime interest. 

Acknowledgments: We thank Luis Anchordoqui, Dario Grasso, Gudlaugur Johannesson, Claudio 
Kopper and Kohta Murase for informative discussions. This work is supported by the U. S. Depart¬ 
ment of Energy under the contract DE-FG-02-95ER40896 and by the National Science Foundation 
under grants OPP-0236449 and PHY-0236449. 

A Galactic Diffuse Gamma-Ray Emission 

The recent measurement of the diffuse Galactic 7 -ray emission by Fermi -LAT [2] covers an energy 
range from 200 MeV to 100 GeV. While we are interested in neutrino energies above 1 TeV, we 
need to have a specific model to perform an extrapolation. Similar to the model used by the Fermi 
Collaboration, we use the GALPR0P code [98-100] to calculate the 7 -ray and neutrino spectra as well 
as their morphology. 

As a starting point, we first reproduce the 7 r° background of the HI gas part in Ref. [2]. Our 
corresponding results are shown in Fig. 11. Note that for the numerical calculations in this paper 
we obtain the total Galactic diffuse neutrino flux from both the HI and molecular gas compo¬ 
nents. Fig. 11 shows the spatial distributions of 7 -ray flux for four different GALPR0P models: 
■ s S’ z 4 i? 20 T 150 c '5, s L z 6 R 20 T ao c 5, s T^10' R 30 t 150 <: 2 and s O z 8 R 30 T co c 2. We use the same short 
hand notation 5 X z z R Rf t Tg c as Fermi -LAT [2] where X is the first letter of the CR source distri¬ 
bution (S: SNR [65], L: Lorimer [ 66 ], Y: Yusifov [101], O: OB stars [102]); Zh and are the vertical 
and radial size of the diffusion region given in units of kpc; T$ is the spin temperature in units of 
Kelvin; c is the E(B — V ) magnitude cut accounting for high-extinction regions in the determination 
of the gas-to-dust ratio. As one can see from Fig. 11, our simulated results can match those from 
Fermi -LAT very well. 

We use GALPR0Pv54_r2504 [103] to model the propagation of the CRs in the Galaxy and calculate 
the photon and neutrino diffused emission from the interactions between the CRs and the interstellar 
medium. For the neutrino flux, we use the energy spectrum and the production cross section 
calculated in Ref. [104]. To be consistent and for the 7 -ray flux calculation, we also use the photon 
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Figure 11: Left panel: a comparison of diffuse 7 -ray distributions in Galactic latitude in the 
inner part of the galaxy with —30° ^ l ^ 30° and 200 MeV ^ E 7 ^ 1.6 GeV for four different 
GALPROP models: s S z 4 R 20 T 150 c 5, s L z 6 R 20 T ao c 5, S Y Z 10 fl 30 T 150 c 2 and s O z 8 R 30 T oo c 2. The 
current black and solid line is from our simulation. The red and dashed line is the HI gas part from 
Fermi -LAT (the top panel of Fig. 20 in [2]). The energy range is from 200 MeV to 1.6 GeV. Right 
panel: a comparison of diffuse 7 -ray distributions in longitudinal latitude for —5° ^ b ^ 5° and 
200 MeV ^ ^ 1.6 GeV. 


energy spectrum and the production cross section in Ref. [104], Our study requires the calculation 
of Galactic CRs above 100 TeV, slightly beyond the recommended energy range of the GALPROP 
code. However, since the morphology of secondary neutrino and photon fluxes only depend on the 
gas target distribution and since we account for uncertainties of the CR spectrum in the knee region 
the additional systematic uncertainties from this extrapolation is not expected to be dominant. 

B Neutrino and Photon Spectra from Dark Matter Decay 

For decaying dark matter models, several models of the energy spectra can provide a reasonable fit to 
the IceCube observed neutrino spectrum [21], As an example, we choose two representative models 
of a DM fermion x and boson X with coupling xHL and X(HL) 2 , respectively, before electroweak 
symmetry breaking. The former operator has the dark matter decay as % —> h + u, while the latter 
operator has the main decay mode of X —> 2h + 2v (see also Ref. [105] for dark matter four-body 
decay). We show a comparison of photon and neutrino spectra for the two representative dark 
matter models in Fig. 12 with the energy bin size to be 10% of the corresponding energy. Compared 
to the observed IceCube neutrino flux, the associated 7 -ray flux can be 10% or above for a wide 
range of dark matter models. For more dark matter models, we show the ratios of photon over 
neutrino fluxes for different dark matter decay modes in Fig. 13. 
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Figure 12: A comparison of the photon and neutrino spectra for two representative decaying dark 
matter models: x —> /i + z/ (left panel) and X —> 2 h + 2v (right panel). The bin size is chosen to be 
10 % of the corresponding energy. 



Figure 13: The ratio of photon over neutrino fluxes for different dark matter decay modes. 
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One could have a neutrino-philic model with suppressed 7 -ray fluxes, although it may not provide 
a good fit to the IceCube observed spectrum. One example is to have the dark matter particle to 
be the electric-neutral component of, X, 3i under SU(2)l x U(1)y- The coupling like L c l XLl can 
induce a two-body decay of the dark matter particle X° to two neutrinos [19]. The suppressed 
three-body decay like X° —> u + e + + W~ can generate a small 7 -ray flux. 

C Background and Signal Sampling 

The background samples of the high-energy starting event (HESE) analysis are produced in the 
following way. The expected zenith distribution of isotropic signal events and atmospheric back¬ 
grounds above a deposited energy of 60 TeV is shown in Fig. 5 of the supplementary material of the 
official IceCube publication [5]. In each of the ten declination bands (equivalent to zenith bands, 
since 9 zen = 7 r /2 + 6) we generate the expected number of events /r for conventional atmospheric 
neutrinos (“atmo z/’), atmospheric muons (“atmo //”) and the best-fit isotropic astrophysical signal 
(“astro”) following an E~ 2 3 spectrum. The zenith angle 9 of each event is uniformly sampled over 
the cos/9 range of the particular bin and the right ascension angle is randomly chosen in [ 0 , 2n). 

For each of the three subsets we generate track (o) and cascade (o) events with a track fraction 
x 0 = 0.9 for atmospheric muons, x 0 = 0.71 for conventional atmospheric neutrinos and x 0 = 0.19 
for the astrophysical signal following the values of Table IV of Ref. [5]. The fraction of cascades 
is denoted by x Q = 1 — 27 • The median angular error 050 % of each track and cascade is randomly 
selected from the actual values reconstructed for the 4 tracks and 16 cascades in the data set 
above 60 TeV (Tab. IV in Ref. [5]). The background events are then re-sampled from their original 
distribution following a von-Mises-distribution. The angular distance a of events from the true 
position is sampled as cos a = 1 + cr 2 ln(l — ,t (1 — exp(— 2 /cr 2 )), where x is a random number in [ 0 , 1 ) 
and the azimuthal direction of the scatter is a random angle in [0, 2n). The angular uncertainty a is 
related to the median provided in Table I in Ref. [5] via cos(* 50 % = 1 + cr 2 ln(l — 0.5(1 —exp(—2/u 2 )). 

The signal samples are generated under the assumption that a fraction 1 — /j so of the best-fit 
astrophysical contribution is due to a Galactic contribution. The first step of the sample generation 
is the same as for the background accept that the isotropic astrophysical expectation is re-scaled 
in each zenith bin to a fraction /i so . For the Galactic distribution we sample events from the 
Galactic signal map w s i gna i which is re-weighted according to the zenith distribution of an isotropic 
component in Fig. 5 of Ref. [5] and then normalized to the expected Galactic event fraction 1 — /i so . 
The distribution between tracks and cascades and the res-sampling of positions in the sky via the 
reconstruction uncertainty is the same as in the case of the isotropic component. 

For the classical v ll: search we are looking for events below zenith angles of 85° to avoid large 
backgrounds from atmospheric muons. The only background that we consider in this case is atmo¬ 
spheric As a benchmark value we expect 10 isotropic astrophysical events in one year with a 
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muon energy larger than about 100 TeV. This is comparable with the 60 TeV energy cut used in the 
HESE analysis. During the same period we expect 10 conventional atmospheric Un events. Signal 
and background events are exponentially suppressed at large zenith angles due to Earth absorp¬ 
tion [106]. Here we account for this absorption via the charged and neutral current interaction of 
100 TeV neutrinos as they traverse the Earth. We use the preliminary Earth density model provided 
in Ref. [106]. The conventional atmospheric v. x is enhanced towards the horizon due to an increased 
atmospheric column depth. We use the parametrization of Ref. [107] to account for this. 

As in the case of the HESE event sampling we generate Galactic signal events by a rescaling of 
the isotropic signal via the fraction f- lso and the corresponding Galactic emission accounting for an 
event fraction 1 — /cal sampled from the corresponding Galactic emission template. Note, that in 
this study we do not have the actual IceCube event distribution for three years of observation and 
we can only use these simulations to study the sensitivity. 

D Likelihood Function 

For the HESE 3yr search the expectation values for tracks (o) and cascades (o) are defined in the 
following way: 

<f 0 (^) = 1 ^^x^y^°(Q)w signa m, (11) 

' A/ Gal 7 

= /iso<f 0 r > astro (H) + > atmo "(«) + xfp V tm ° "(«), (12) 

where ,t^ 0 corresponds to the fraction of tracks and cascades for the three separate components 
and fj, are the expectation values extracted from Fig.5 and Tab. IV in Ref. [5[. The factor 
renormalizes the Galactic flux such that the total number of expected events is conserved, 

A/cai = JdO/<r(0)^ signa i(Q) (Jdfi^(fi)) 1 • (13) 

Numerically, we find that A/Qai is in the range 1.02 (Galactic diffuse) to 1.09 (DM decay) except for 
the Fermi Bubble template where Mq & \ — 1.24. 

The three year expectation values of the isotropic astrophysical contribution (Mo/!™), the atmo¬ 
spheric neutrino ( v ) and the atmospheric muon >l ) are known from Fig. 3 and Tab. IV 

of Ref. [5[. Now, the angular uncertainty of the events can be accounted for by replacing the 
expectation values for signal and background by the weighted values 

/h = Jdf2pdfj(J7) /x(H). (14) 

The fraction X{ indicates the contributions into tracks and cascades. For normalized pdf’s this 
implies that /hot = //tot and the exponential terms in the likelihood drop out of the likelihood ratio 
in Eq. (9). 
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The case of the classical search follows the same prescription. However, due to the good 
angular resolution of tracks we choose an initial binning of the skymaps such that each pixel has an 
approximate diameter of 1°, corresponding to about twice the angular resolution of IceCube. In this 
case we do not need to resample the angular resolution from the true position. Also, we do not have 
to account for different event topologies in the simulation, i.e. x 0 = 1. The background and signal 
expectation are both reduced towards high zenith angles due to Earth absorption via charged and 
neutral current interactions. The conventional atmospheric flux is enhanced towards the horizon 
due to the scaling with the atmospheric depth. We already commented about this approximation in 
the previous section. The total expected astrophysical z/ ;i event rate is normalized to 10 events per 
year and the total expected conventional atmospheric u iL event rate is also normalized to 10 events 
per year. This approximates the expected integrated event rate above muon energies of 100 TeV [85]. 
Due to the limited field of view with 9 zen > 85° the normalization factor (13) is smaller than unity: 
0.79 (PWN), 0.81 (DM decay), and 0.86 (SNR & Galactic diffuse). 


E Comparison of Galactic and Extragalactic Emission 


Any mechanism that provides a strong anisotropic Galactic neutrino emission is also expected to 
contribute via an isotropic emission from similar, but distant galaxies. If Q V (E) is the total spectral 
emission rate of our Galaxy then the contribution of extragalactic source can be written as the 
red-shift integral [108] 




I 


dz 

H(z) 


p{z) Q v [( 1 + z)E v ] , 


(15) 


where H(z)/Hq = a/Qa + H m (l + z ) 3 is the red-shift dependent Hubble constant in the matter- 
dominated era with Hq = 67.27 + 0.66 kms -1 Mpc -1 and O m = 1 — = 0.3156 + 0.0091 [109]. 

The distribution of sources is taken to be its averaged value p( 0) = po ~ (10 -3 — 10 -2 ) Mpc -3 . 
Comparing to the average diffuse emission from the Galaxy can be expressed via Eq. (6) and 
approximating the neutrino flux as power-law with index T, we have the relation 

cp{z) (1 + z)~ r cpoCz „ in _ 3 . f Po \ /(20kpc) -2 
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with source evolution factor = 0(1) (£ 2 ^ 0.5 for no evolution and T ~ 2). 

Note that the parameter J is defined as ((r — rQ) -2 )/(47r), averaged over the source distribution 
in the Galaxy. The SNR and PWN distributions in Eq. (7) are essentially 2-dimensional and have 
a maximal contribution at r max = olRq/^. For r max < Rq the size of the J-factor is then expected 
to be of the order of 47 tJ ~ 0((Rq — r^ ax ) -1 ]. However, this does not take into account the 
radial width of the distribution. Numerically, we find 47 tJ ~ (4.6 kpc) -2 for both, SNR and PWN 
distributions. The local Galactic emission is then typically much stronger than the corresponding 
extragalactic contribution. 
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On the other hand the dark matter distribution in Eq. (8) is 3-dimensional with a characteristic 
extension r s » Rq. The dark matter profile is extended to a wider range compared to the SNR 
or PWN sources, and therefore has a larger value of J -1 / 2 . The characteristic length scale of the 
J-factor is of the order of 4irJ ~ 0(rJ 2 /S). This agrees well with the numerical value 4-7 tJ ~ 
(34kpc) -2 . For a flatter dark matter profile, i.e. larger r s , we also anticipate a larger value of J^ 1 / 2 
or a smaller value of J and more extragalactic contributions. In the case of dark matter decay the 
local density in Eq. (16) is given by po — ^cdmPcr/^DM al — 1CT 2 Mpc -3 [109]. Therefore, up to 
uncertainties of the Galaxy density and the dark matter distribution, the isotropic emission from 
extragalactic dark matter is expected to be at a similar level as the Galactic emission [20,25]. 
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